function analyzeDistributionsForTherma(pattern,workspace)

eval(['load ' workspace]);

variables = who(pattern);

thermRate = 0:1/(length(variables)-1):1;

for i=1:length(variables)
    data = eval([char(variables(i))]);
    mean(i) = mean(data);
    sigma(i) = std(data);
    [H(i),pValue(i)] = jbtest(data);
    
end

subplot(2,2,1);
plot(thermRate,mean,'.k');
ylabel('mean');

subplot(2,2,2);
plot(thermRate,sigma,'+r');
ylabel('std');

subplot(2,2,3);
plot(thermRate,pValue,'.k');
ylabel('p-value');

subplot(2,2,4);
plot(thermRate,H,'.k');
ylabel('test de Hipotesis');